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Abstract: Supersymmetric grand unified models based on the gauge group SO (10) are 
especially attractive in light of recent data on neutrino masses. The simplest SO (10) 
SUSY GUT models predict unification of third generation Yukawa couplings in addition 
to the usual gauge coupling unification. Recent surveys of Yukawa unified SUSY GUT 
models predict an inverted scalar mass hierarchy in the spectrum of sparticle masses if the 
superpotential fi term is positive. In general, such models tend to predict an overabundance 
of dark matter in the universe. We survey several solutions to the dark matter problem in 
Yukawa unified supersymmetric models. One solution- lowering the GUT scale mass value 
of first and second generation scalars- leads to ur and cr squark masses in the 90—120 GeV 
regime, which should be accessible to Fermilab Tevatron experiments. We also examine 
relaxing gaugino mass universality which may solve the relic density problem by having 
neutralino annihilations via the Z or h resonances, or by having a wino-like LSP. 
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Standard Model. 



1. Introduction 



The unification of gauge couplings in the Minimal Supersymmetric Standard Model (MSSM) 
can be construed as indirect evidence for the existence of weak scale supersymmetry. It 
is of course also suggestive of the existence of a grand unified symmetry at energy scales 
Q > Mqut - 2 x 10 16 GeV. Such a SUSY GUT theory may be the "low energy effective 
theory" that can be obtained from some more fundamental superstring theory. Models 
based on the gauge group SU(5) are compelling in that they explain the apparently ad-hoc 
hypercharge assignments of the SM fermions|]l|. However, many SU(5) SUSY GUT mod- 
els as formulated in four dimensions are already excluded by proton decay constraints Q. 
SU(5) SUSY GUT models can also be formulated in five or more dimensions, where com- 
pactification of the extra dimensions leads to a break down in gauge symmetry ||]. These 
models can dispense with the unwieldy large Higgs representations required by four dimen- 
sional models, and can also be constructed to suppress or eliminate proton decay entirely. 

The steady accumulation of data on neutrino mass and neutrino oscillations in the 
past several years has given strong motivation to the possibility that nature is described 
by a SUSY GUT theory based on the gauge group 5O(10)[|, g. In 50(10) theories, 
all 15 of the SM matter fields of a single generation inhabit the 16-dimensional spinorial 
representation of 50(10). The 16th element can be taken to be a SM gauge singlet field 
that functions as the right-hand neutrino. A large Majorana mass term is allowed for this 
field by the gauge symmetry. Dominantly left handed neutrinos thus gain tiny masses 
via the see-saw mechanism, while dominantly right handed neutrinos acquire a mass near 
to the unification scale, and are conveniently swept out-of-sight Q . The neutrino see-saw 
mechanism also lends itself elegantly to a mechanism for baryogenesis via leptogenesisQ. 
For these reasons (and others), the possibility that nature is described by an 50(10) SUSY 
GUT model at scales Q > Mqut deserves serious consideration. 

In the simplest 50(10) SUSY GUT models, the two Higgs doublets of the MSSM 
occupy the same 10 dimensional Higgs multiplet 0(10). The superpotential then contains 
the term 

/ 9 M16M16M10) + ■ ■ ■ (1.1) 

where / is the single Yukawa coupling for the third generation. Thus, the simplest 5O(10) 
SUSY GUT models predict t — b — r Yukawa coupling unification in addition to gauge 
coupling unification. It is possible to calculate the t, b and r Yukawa couplings at Q = 
n^weak, and extrapolate them to Mqut m much the same way one checks the theory 
for gauge coupling unification. Yukawa coupling unification turns out to depend on the 
entire spectrum of sparticle masses and mixings || since these enter into the weak scale 
supersymmetric threshold corrections. Thus, the requirement of t — b — r Yukawa coupling 
unification can be a powerful constraint on the soft SUSY breaking terms of the low energy 
effective theory |]]. 

It is well known from early work that large values of tan (3 ~ 50 are required for 
t — b — t unification. However, in models with a common scalar mass, such as mSUGRA, 
it was found that a breakdown in the mechanism of radiative EWSB occurred which ex- 
cluded tan/3 values which would give rise to Yukawa coupling unification. In Ref. [10|, it 
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was found that models with Yukawa coupling unification with R < f.05 (good to 5%) 
could be found if additional D-term splittings of scalar masses were included. Here, 
R = max(ft, fb, f T )/ mini ft, fb, fr), where all Yukawa couplings are evaluated at the 
GUT scale. The D-term splittings occur naturally when the 50(10) gauge symmetry 
breaks to SU(5), and they are given by[[Tl|] 



m Q = m E = m u = m w + M D> 
m D= m L = m 16 - 3M Di 

m 2 N = m\ G + 5M|), 
,2 



m z Hud = mi T2M z D , (1.2) 

where parameterizes the magnitude of the D-terms. Owing to our ignorance of the 
gauge symmetry breaking mechanism, can be taken as a free parameter, with either 
positive or negative values. \Md\ is expected to be of order the weak scale. Thus, the 
D-term (DT) model is characterized by the following free parameters, 

mi 6 , mio, Mf,, m 1/2 , A , tan/?, sign(/i). (1.3) 

Using the DT model, Yukawa unification good to 5% was found when soft term parameters 
mis and m-io were scanned up to 1.5 TeV, but only for fj, < values [[lC|]. The essential 
quality of the D-term mass splitting is that it gave the value of mn u a head start over mn d 
in running towards negative values, as is required for REWSB. In Ref. p2[ , the DT model 
for \x < was explored in more detail, and the neutralino relic density was explored. It was 
found that a good relic density could be found in the ^4-annihilation funnel p3|, although 
a rather large value of BF(b — > 57) ~ 4 x 10" 4 was usually predicted. Good relic density 



was also found in the stau co-annihilation [14] and hyperbolic branch/focus point (HB/FP) 



region [15, 16, [P?]], but at some cost to the degree of Yukawa coupling unification 1 . 

Later, E821 measurements of the muon anomalous magnetic moment appeared, which 



favored positive values of the superpotential [i parameter. In Ref. [18|, it was found that 
Yukawa coupling unification good to only 30% could be achieved in DT models with fi > 
when 77116 values were scanned up to 2 TeV. The models with the best Yukawa coupling 
unification were found to have soft term relations 

Al ~ 2m 2 w ~ 4m? 6 , (1.4) 

which had also been found by Bagger et al. in the context of radiatively driven inverted 
scalar mass hierarchy (IMH) models ||, ||, || 2 . In Ref. f|], a model with J ust GUT scale 
Higgs mass splittings was also examined (HS), while all other scalars remained universal. 



The parameter space of the HS model is that of Eq. 1.3, but where the D-term splitting is 



1 The HB/FP region is characterized by \fi\ ~ M2 so that the LSP is a mixed higgsino-bino, and has a 
large annihilation rate into vector bosons. 

2 In IMH models, starting with multi-TeV scalars masses at Q — Mgut, for certain choices of boundary 
conditions third generation and Higgs soft masses are driven to sub- TeV values, while first and second 
generation scalars remain heavy. 
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only applied to the last of the relations in Eq. L2. Yukawa coupling unification in the HS 
model was found to be comparable to the DT model case when m\Q values up to 2 TeV 
were scanned. Finally, in Ref. ]22| (Auto et al), soft term values of m\§ up to 20 TeV were 
explored. In this case, Yukawa unified solutions to better than 5% were found for fi > for 
the HS model when very large values of m\§ > 5 — 10 TeV were scanned. The large scalar 
masses that gave rise to Yukawa unification also acted to suppress neutralino annihilation 
in the early universe, so that rather large values of the relic density were found. Models 
with f^/i 2 < 0.2 could be found, but only at the expense of accepting Yukawa coupling 
unification to 20%, rather than 5%. The models found with low relic density generally had 
either a low [i value, or were in the light Higgs annihilation corridor, with 2m^ i ~ m^. 

Similar work has been carried on by Blazek, Dermisek and Raby (BDR). In Ref. 
1 23, 24], the BDR group used a top-down approach to the RG sparticle mass solution to 
find Yukawa unified solutions for [i > 0, where they also noted that in this case the HS 
model worked better than the DT model. In their approach, the third generation fermion 
masses and other electroweak observables were an output of the program, so that starting 
with models with perfect Yukawa coupling unification, they would look for solutions with a 
low x 2 value constructed from the low energy observables. The BDR Yukawa unified solu- 
tions were also characterized by soft term IMH model boundary conditions. The solutions 
differed from those of Ref. |^2| in that they always gave a very low value of Higgs mass 



rriA and also small fj, parameter, indicative of a mixed higgsino-bino LSP. In Ref. [25], the 
neutralino relic density was examined for the BDR solutions. Their low \x and tha values 
generally led to very low values of Q^h? unless mi/ 2 was small enough compared to [i that 
the LSP was in the mixed higgsino-bino region. 3 

In this paper, we give further scrutiny to the Auto et al. Yukawa unified solutions, 
especially with respect to the relic density of neutralinos expected from these models. 
Our goal is to find regions of parameter space where a high degree of Yukawa coupling 
unification is preserved, but also where a cold dark matter relic density in accord with 
WMAP measurements p7| can be found. 4 To do so, we expand the parameter space of 
the SO( 10) model. In Sec. ||, we consider relaxing generational universality of the matter 
scalars. The expansion of parameter space in this fashion is consistent with the SO( 10) 
gauge symmetry, which only preserves universality within a generation (when SO( 10) is 
unbroken). By maintaining TeV scale third generation masses, but reducing first and 
second generation masses, light first and second generation scalars are produced which 
enhance the neutralino annihilation rate. The curious effect in the HS model is that the 
split Higgs mass contribution to RG running causes the right up and charm squarks to be 
by far the lightest of the scalars, with m% R and rric R ~ 100 GeV needed to achieve the 
required relic density. We show that such light squarks are just at the limit of observability 
for the Fermilab Tevatron collider. Our analysis motivates an experimental search for 
just two species of light squarks, while all other scalars are in the TeV range, and where 
m-g ~ 400 GeV. 

In Sec. ||, instead we relax the condition of gaugino mass universality. The gauge kinetic 
function of supergravity theories need not adopt its minimal form in order to respect the 
SO (10) gauge symmetry. We examin e several models found in Ref. p^| , and vary the U(l) 



3 We note that the HS model without Yukawa coupling unification has recently been examined by Ellis 
et al, Ref. Q. 

4 The WMAP collaboration has determined Q C DAih 2 = 0.1126±omli at 2a level. 
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gaugino mass M\. This allows the bino-like LSP to sit atop the h or Z resonance, which 
in some cases is enough to reduce the relic density to the ^j^/i 2 ~ 0.1 level. Also, by 
increasing substantially the value of Mi, ultimately the lightest neutralino becomes wino- 
like, and a valid relic density is also achieved. Variation of the gaugino masses preserves 
the Yukawa coupling unification in all the models we examined. We also examine the case 
of lowering the SU(2) gaugino mass M 2 to obtain a wino-like LSP. In this case, the lightest 
chargino usually drops below limits from LEP2. In Sec. we present a summary and 
conclusions. We note here additional studies which have been performed, where just b — r 
unification was examined) 25, 29] or where Yukawa quasi-unification was examined 



2. Scalar mass non-universality and light squarks 

Our goal in this section is to explore whether an expanded parameter space allowing non- 
universality of generations can provide new avenues towards solving the dark matter prob- 
lem in Yukawa unified supersymmetric models. Thus, we expand the parameter space of 
the HS model to now include 

mi 6 (l), mi 6 (3), mio, Mp, m 1/2 , A , tan (3 and sign (//), (2.1) 

where the matter scalar soft mass is now broken up into mig(3) for the third generation, 
and mi6(l) for the first generation. We will also assume mi6(2) = miQ(l) to reduce 
parameter space freedom and to quell potential contributions to flavor changing processes 
like K^ — Kg mass difference and BF(u — > ej). Allowing third generation non-universality 
can also contribute to the Amg = m® B — m B mass difference, the rate BF{t — > cy) or 
BF(t — > U7), and BF(b — ► S7). However, these latter constraints are all rather mild, and 



in general can allow for significant deviations in generational non-universality [31, 32 



Our goal is to maintain the successful prediction of Yukawa coupling unification found 



in Ref. [22] with heavy scalars and an inverted scalar mass hierarchy, but at the same time 
reduce some of the scalar particle masses so that there is not too severe a suppression 
of neutralino annihilation in the early universe which leads to too large a relic density. 
The IMH mechanism of Bagger et al. works if the mass relation between third generation 
scalars, Higgses and A terms is preserved: A^ ~ 2m 2 ~ 4m 2 6 (3). Thus, a possible way 
forward is to decrease the value mi6(l) = mi6(2) so far that a valid relic density is attained. 
It turns out this procedure is possible, but that it leads to a qualitatively new spectrum 
of superparticle masses wherein two of the squarks- ur and cr- become quite light, in the 
90-120 GeV regime, while the other scalar masses remain quite heavy. 

A sample spectrum is shown in Table [l], which uses Isajet 7.69 to generate the spectrum|34 
As can be seen, a high degree of Yukawa coupling unification is maintained, while simul- 
taneously achieving an acceptable dark matter relic density ^j^ 2 = 0.09. The low relic 
density is obtained because the light ur and cr squarks allow Z±Z± — > uu, cc annihilations 
to proceed at a high rate due to unsuppressed t-channel squark exchange. 

Why are some of the squark masses quite light, while most other scalars are heavy? 
The answer comes from the 1-loop RGEs for the SUSY soft breaking terms. For first 
generation matter scalars, they have the form j3f| 
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where 



+ Tr [ 



Hi/ 



m| 



2m^ + mf, + 



In the mSUGRA model, the term involving 5 
doesn't contribute to RG running at all, since uni- 
versality implies 5 = 0. However, in the HS model, 
while the trace over mass matrices remains zero, the 
contribution from Higgs splitting is large, especially 
for our case where the soft Higgs masses can be ~ 10 
TeV. Radiative EWSB in Yukawa unified models re- 
quires fn 2 H d > m ?/„! so that in fact S < 0. Then the 
soft masses with negative co-efficients on the S terms 
will be highly suppressed in their RG running, while 
terms with positive S coefficient will be pushed to 
higher mass values. If the S term dominates the RG 
running, then indeed we expect a large splitting due 
to RG evolution in all the first and second genera- 
tion soft breaking masses, with the splitting being 
dictated by the soft mass U(l)y quantum numbers. 
The splitting amongst third generation soft masses 
will be much less, since their RGEs include addi- 
tional Yukawa coupling times soft mass terms, and 
the third generation soft terms are also large. 

In Fig. [l]a., we show the third generation and 
Higgs soft term running, for the same case as in Ta- 
ble |]. The soft term suppression due to Yukawa 
coupling contributions is apparent. In Fig. [l|&., we 
show first generation soft term running. In this case, 
the large splitting due to the dominant S term is evident, and causes the U\ mass to run 
to very low values, which gives rise to the light ur and cr squarks. In addition, the L\ 
soft mass also is suppressed, though not as much as U\, so that rather light left selectrons 
and smuons and their associated sneutrinos are generated. 



parameter 


value (GeV) 


, %e(3) 


7826.5 




1202.0 


mio 


9652.3 


8$) 


(2894.9) 2 


mi/2 


80.0 




-16626.0 


tan (3 


51.165 




0.563 


fb(M GUT ) 


0.549 


U(Mgut) 


0.563 


R 


1.02 


#!) 


3075.4 


m- g 


380.1 




1172.2 




107.0 




1175.0 


^ 


1343.7 


m iL 


693.4 


m iR 


1734.2 


m H 


672.7 


m h 


1651.1 


mr 


2330.5 




3049.2 


m m 


135.1 




57.3 


m z 2 


134.9 


m A 


2505.7 


m h 


134.4 


n z 1 h2 


0.093 


BF(b -» 87) 


3.96 x 10~ 4 




15.7 x 10~ 10 



Table 1: Masses and parameters 
in GeV units for a Yukawa uni- 
fied HS model with generational non- 
universality. The spectrum is obtained 
using ISAJET v7.69. 
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In Fig. ^, we show the mie(l) vs. mj/2 plane for fixed mie(3) = 7826.5 GeV, and 
other parameters as given in Table |. The blue region at low m^2 is excluded by the 
LEP2 bound on chargino mass > 103.5 GeV. The black region on the left is excluded 
because the rny 1 squared mass is driven tachyonic, resulting in a color symmetry violating 
ground state. The green region denotes sparticle mass spectrum solution with relic density 
0.094 <VL^h 2 < 0.129, within the WMAP favored regime, while the yellow region denotes 
even lower relic density values, wherein the CDM in the universe might be a mixture 
of neutralinos plus some other species. The red region has % h 2 < 0.5. It may be 
allowed if QCD corrections to Z\Z\ — » uu, cc are large and positive 5 . The remaining 
unshaded regions give sparticle spectra solutions with h 2 > 0.5, and would be excluded 
by WMAP. We also show several contours of BF(b — ► sj), Aa'^ USY and R, the max/min 
ratio of GUT scale Yukawa couplings. We see that all constraints on the sparticle mass 
spectrum are within allowable limits for the shaded region to the right of the excluded 
region. 

We have also performed a x 2 analysis in the parameter space listed in Eq. 2.1 and 
scanned over the following range of seven parameters: 

< mie(l) < 10 TeV, < m 16 (3) < 10 TeV, < m 10 < 15 TeV, 

< m 1/2 < 0.5 TeU, < M D < 5 TeV, -30 < A < 30 TeV, 40<tan^<60. (2.8) 

For every point in parameter space, we calculate the x 2 quantity formed from Aa^, BF(b — ► 
sj) and Qgji 2 . We present plots with a color coded representation of the x 2 value, where 
red shading corresponds to high x 2 , green shading corresponds to low x 2 and yellow shading 
has intermediate values. The numerical correspondence is listed in Fig. [||. As a result, 
we have found that in order to have Yukawa unification below ~ 5% level and low x 2 , one 
should be in the special, strongly correlated region of the parameter space. For further 
analysis we have fixed some parameters at the values favored by the x 2 analysis: 

m 16 (3) = 7750 GeV, A = -2.05 m 16 (3), (2.9) 

and varied the others within the x 2 fit preferred range: 

49.8 < tan [3 < 51.6, 75 GeV < m 1/2 < 120 GeV, 

0.14 < mH * ~ mg " < 0.23, 1.225 < < 1.27. (2.10) 

mi 6 (3) 2 mi 6 (3) 



Finally, for each point generated in the range specified by Eq. |2.10 the value of \ 2 wa s 
minimized by varying mi6(l). 

Those solutions and respective correlations among the model parameters are illustrated 
in Fig. |3] where we show the results of scanning over the complete parameter space listed in 
Eq. |2.8| . The color coding of the points corresponds to a x 2 value computed from the value 
of Jl^/i 2 , Aa^ USY and BF(b — > sj). In constructing the x 2 quantity, we adopt the limits 

5 A subset of QCD corrections for ZiZi — > qq has been calculated by Drees et al, and are found to be 
small H. 
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0.1126 ±0.00805 (we use the upper limit only), Aa^ USY 



(27 ±10) x 10~ 10 
measurement compared to SM predictions |57]]), 



that h 2 

(as required by the recent E821 (g — 2) M 
and BF(b -> 57) = (3.25 ± 0.54) x 10" 4 |8]. 

In Fig. ||, we show the results of scanning over the narrowed parameter space defined 
by Eq. 2.1C in the R vs. ^^Ji 2 plane with the same color coding as in Fig. |3| We see that a 
large range of points are generated with low x 2 value, and that many of these have Yukawa 
coupling unification below the 5% level, and also with a WMAP allowed relic density. 




Q (GeV) 

Figure 1: Evolution of a.) third generation and Higgs scalar soft masses and b.) first generation 
soft masses for the case of m 16 (3) = 7826.5 GeV, m w = 9652.3 GeV, M D = 2894.9 GeV, m w {\) = 
1202.3 GeV, m 1/2 = 80 GeV, A a = -16626 GeV, \i > 0, tan/3 = 51.165 and m t = 180 GeV. For 
Higgs soft masses, we actually plot sign(rrvjj) ■ ^/|m^-|. 



In Fig. [|, we show the same points from the parameter space scan, but this time in 
the m,i„ vs. R plane. In this case, we see that essentially all of the low x 2 points have 



m 



values in the 90-120 GeV range. For models with m„ B > 120 GeV, the value of ^-^^h 2 



exceeds the WMAP bounds, and the models become excluded. The allowed models shown 
all have R values in the 1.02-1.21 range. 
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HS SC(10) model: m ie (3}=7826.5GeV. m i0 =9652.3GeV, M^n\ 16 (3)=0.3S988 




40 



H/fV 0.094 
0.094< Llh*< 0.129 
0.129< &tf< 0.5 



LEP2 

■ 



1.16 1.17 1.18 1.19 1.2 1.21 1.22 1.23 1.24 

Figure 2: Allowed parameter space of YukCl?af $firftli fvrp'Ss|fjiimetric HS model with generational 
non-universality. We show the mie(l) vs. plane for mie(3) — 7830 GeV, mio = 9650 GeV, 

M D /m 16 {3) = 0.37, A Q /m 16 {3) = -2.1, fi > 0, tan/3 = 51 and m t = 180 GeV. The black 
shaded region gives tachyonic particles, while the blue region is excluded by LEP2 chargino search 
experiments. The yellow and green regions are allowed by the WMAP determination of fi^/i 2 . We 
also show contours of R, the measure of Yukawa unification at Mqut- 



2.1 Phenomenological consequences of two light squarks 

Given that the Yukawa unified model with split generations predicts two light squarks of 
mass 90-120 GeV in order to obtain a relic density in accord with WMAP, it is natural to 
examine the phenomenological consequences of such light squarks. We first turn to data 
from LEP2. In Ref. the ALEPH collaboration reports mass limits on squarks in 

the MSSM. Using 410 pb _1 of data ranging up to CM energies of 201.6 GeV, they report 
rriq > 97 GeV at 95% CL assuming q — > qZ\ with a mass gap rriq — > 6 GeV, and 
assuming 10 species of degenerate squarks. In our case, with just two nearly degenerate 
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squarks, the mass limits will be even lower 6 . Furthermore, the OPAL collaboration^] has 
searched for t\t\ pairs assuming t\ — ► cZ\. These limits exclude can be applied to our case 
of crcr production, and exclude < 90 GeV for almost all values of < ra^ R — m c . 

There are also limits on squark pair production from Fermilab Tevatron experiments. 
In Ref. [41 1, the CDF collaboration reports a mass limit of rriq ~ 250 GeV when m~ g ~ 400 



GeV, based on 84 pb~ x of data taken in Run 1 with ^fs = 1.8 TeV. This analysis assumes 
eight degenerate squarks when MSSM parameters are used, and can potentially rule out 
the case presented here with two species of nearly degenerate squarks. Thus, we examine 
the relative signal produced in the Yukawa unified case with two light squarks with mass 
~ 90 — 120 GeV. A similar analysis by the DO collaboration [E^] also excludes rriq < 250 
GeV (with 10 degenerate squarks) by searching for jets+ $t events, using 79 pb _1 of Run 
1 data. This analysis assumes the mSUGRA model, however, which doesn't give rise to 
our case where rriq <C rrig. 

We generate all sparticle pair production reactions using the event generator Isajet 
7.69, using the model parameters as given in Table |], except allowing m\§{l) to float as a 
free parameter. The mass value varies nearly linearly with mi6(l), so this allows us 
to vary the squark masses, while keeping other parameters fixed. We use the ISAJET toy 
detector CALSIM with calorimetry covering the regions —4 < r\ < 4 with cell size Ar] x 
A(f> = 0.1 x 0.262. Electromagnetic energy resolution is given by AE em /E em = 0.15/ \JE em , 
while hadronic resolution is given by AE^/E^ = 0.7/\/£a- Jets are identified using the 
ISAJET jet finding algorithm GET JET using a fixed cone size of AR = ^ Ar] 2 + A(f> 2 = 
0.7. Clusters with E? > 15 GeV and \rj(jet)\ < 2.5 are labeled as jets. Muons and electrons 
are classified as isolated if they have > 5 GeV, \r]i\ < 2.5, and the visible activity within 
a cone of R = 0.4 about the lepton direction is less than 2 GeV. 

The CDF analysis requires the following cuts: 

• N(jets) > 3, 

• at least one jet has \r]\ < 1.1, 

• no isolated leptons, 

• the second and third highest Et jets are fiducial, 

• Ri = + (tt - 6<f>() > 0.75 rad, R 2 = y/6</>( + (vr - 5^ 2 ) > 0.5 rad (Sfa = 

pleading jet ~ 4fy T \ an d = | ^second jet ~ <ttyr\), 

• Et > 70 GeV, 

• H T =E T + E T (jet 2) + E T (jet 3) > 150 GeV, 

which we also implement into our analysis. In Fig. |(| we plot the signal cross section after 
cuts. The dashed curve is the result including just qq + gq + gg production (as CDF does), 



6 If we make the simplifying assumption that all of the mass limit comes from running at the highest 
energy of y/s = 201.6 GeV, then the 97 GeV mass limit for 10 degenerate squarks would be reduced to 89 
GeV if only urUr and crcr are produced. 
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while the solid curve includes all sources of > 3 jet events, including especially W^W^ 
and W^Z2 production. The points and error bars correspond to our actual signal rate 
computation, along with statistical error. We also show the 95% CL limits as computed 
by adding maximally allowed signal cross section with SM error either linearly (BG2) or 
in quadrature (BG1). We see that, depending on how the signal and background levels 
are interpreted, squark masses in the 100 GeV range may or may not be allowed by Run 
1 analyses. For the case shown, niu R ~ 100 — 120 GeV would be excluded. If instead we 
take a case with larger m 1/2 = 99.5 GeV, m w (l) = 1297.1 GeV, m w (3) = 8767.2 GeV, 
mio/mi 6 (3) = 1-238, M D /m m (3) = 0.356, tan/3 = 50.79 and A /m 16 {3) = -2.13, then 
we generate a model with R = 1.02, £1% h 2 = 0.13, m„ fl = 117 GeV but with heavier 
charginos and neutralinos, with <r(3 — jets+ $t) = 203 fb, which is below both BG1 and 
BG2 levels. In Fig. ^, we show model points in the vs. rriu R plane. In this case, we 
see the mass gap rriu R — ranging from 25-50 GeV. (The bands are due to sampling of 
mi/2 m steps of 5 GeV, while other parameters are varied continuously.) Solutions with a 
large mass gap are more likely to give large observable rates at Tevatron searches, while 
solutions with small mass gaps will be more difficult to see, since there is less visible energy 
released in the squark decays. 

An overall assessment seems to be that the two light squark scenario from Yukawa 
unified models with generational non-universality is at the edge of exclusion from previous 
Tevatron searches using Run 1 data. However, given that Run 2 has already amassed 
~ 400 pb _1 of data, a new dedicated analysis by CDF and/or DO of > 2— jets + $t events 
arising from squark pair production will likely either rule in or rule out this scenario with 
two light squarks with mass in the 90-120 GeV range. 



We mention here that there also exist Tevatron searches for pp — > tit±X production! 43], 
where it is assumed that i\ — > cZ\ 100% of the time. These analyses by CDF [44] and DO [45] 
require a tagged charm jet, so would not apply to urur production. They do apply to 
CRtcR production. However, the best limit (by CDF) requires < 50 GeV, which is too 
light to apply to the scenario given here: see Fig. 0. However, a new search for cc+ j^T 
events by CDF and DO using Run 2 data should be able to probe the parameter space of 
the scenario listed in this paper. 

Another consequence of light squarks is that they can mediate neutralino-nucleon scat- 
tering at observable rates via squark exchange in s and u-channel quark-neutralino sub- 
process diagrams. In Fig. ||, we plot various 50(10) HS models with Yukawa unification 
to better than 20%, in the log 10 a(Zip) vs. plane. We use the neutralino-nucleon 

scattering rate code developed in Ref. |4(| . The colors of the dots correspond to the \ 2 
value as presented in Fig. |I[ We see that the models with low relic density falling within 
the WMAP acceptable range also have rather large values of neutralino-nucleus scattering 
cross section. In fact, the purple (blue) dots are already excluded by Edelweiss |4?J (CDMS 
II [fll|) searches. CDMS II is expected to ultimately probe neutralino-proton scalar cross 
sections as low as 10~ 8 pb. This essentially covers the complete range of models with 
acceptable Q^^ 2 as obtained from Yukawa unified models with light uji and cr squarks. 



- 10 - 



3. Gaugino mass non-universality 



Aside from generational non-universality in 50(10) models, it is also possible to have non- 
universality of gaugino masses. In a general non-renormalizable SUSY gauge theory, the 
gaugino masses arise from the term 



C 3 -i J d 2 9 L f A B(S)W%W B , (3.1) 



where we have introduced the gauge kinetic function (GKF) /ab(S) which is an analytic 
function of chiral superfields Si, and where Wa are gauge superfields. In supergravity 
theories, the GKF can be expanded as 

fAB = 5 A B + j^ + ---, (3.2) 

where, to maintain gauge invariance, 4>ab transforms as the symmetric product of two 
adjoints. If the auxiliary field F^ab of 4>ab obtains a SUSY breaking (and possibly gauge 
symmetry breaking) vev, then gaugino mass terms are induced: 

C 3 {I 0-^a\b + ■■■■ (3.3) 



In the case of SU(5) GUT theories |5C], the product of two adjoints can be decomposed as 



(24 x 24) sym = 1 + 24 + 75 + 200, (3.4) 

where only the singlet represention leads to universal gaugino masses. The higher repre- 
sentations lead to relations amongst the GUT scale gaugino masses which are dictated by 
group theory, and are listed in Ref. [5C]. In the general case, the field ipAB could trans- 
form as a linear combination of the above representations, in which case the three MSSM 
gaugino masses at Q = Mgut would be essentially free parameters. In the case of 50(10) 



GUTS 13, H], we have 

(45 x 45) sym = 1 + 54 + 210 + 770. (3.5) 



For our purposes here, we will not adopt any specific GUT symmetry breaking pattern 
or set of gaugino masses, but instead merely enlarge the parameter space by allowing each 
MSSM gaugino mass as a free parameter. Given this freedom, it is then possible to take any 
of the Yukawa unified solutions generated by Auto et al. p2|| and vary the gaugino masses to 
see if an acceptable relic density solution is found. The only problem of concern is whether 
doing so will destroy the Yukawa coupling unification, or lead to models which are violated 
by constraints such as those imposed by LEP2 searches for sparticles. We note here that 
an alternative approach is to allow restricted, or even complete freedom of gaugino masses, 
and to explore whether Yukawa coupling unification is possible. That approach has been 
followed in the context of 5-d gaugino-mediated SUSY breaking models fl53"f , where all scalar 
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masses are set to zero at M c ~ McuTi we refer to the work of Balazs and Dermisek for 
more details [54]. 



We have examined the positive fi Yukawa unified solutions presented in Table 1 of [22] 
with regard to variation of the GUT scale gaugino masses M\ and M.2- In the case of 
decreasing Mg, ultimately the lightest neutralino makes a transition from being bino-like 
to being wino-like. Wino-like LSPs have a large annihilation and co- annihilation rate into 
SM particles, and give a relic density almost always much lower than the WMAP bound] 



However, in lowering M2 to the point of achieving a valid Q?,Ji 2 , we find in these five cases 
that the chargino mass always falls below the LEP2 limits, so that the model becomes 
excluded by LEP2 rather than WMAP. The other case, that of varying the U(l)y gaugino 
mass Mi, leads to models that can satisfy the WMAP bound while being consistent with 
other constraints including those from LEP2. 

In Fig. pi, we show the relic density for four of the points from Table 1 of p2| versus 
variation in GUT scale gaugino mass M\. In all cases there is a similar structure: for 
generic values of Mi, the relic density lies beyond the WMAP limits. However, there exists 
a double dip structure where the relic density may or may not fall below the WMAP limit, 
and then at high Mi, the relic density drops to very low values which would imply non- 
neutralino cold dark matter. The high Mi behavior occurs because as Mi is increased, 
ultimately Mi > M2, and the LSP becomes wino-like, although in this case, since M2 is 
fixed, the chargino mass stays fixed beyond the LEP2 mass limits. At lower Mi values, 
the double dip structure occurs when the value of 2m^ i drops below the h or the Z 
resonance. In these cases, the neutralinos can efficiently annihilate at high rates through 
these resonances in the early universe, leading to lower relic density values. In some of the 
cases, the diminution is sufficient to drop below WMAP limits, while in others (those with 
very heavy scalars, and hence with a large /j parameter), the diminution is insufficient to 
save the model. In each case, we have found that the Yukawa coupling unification changes 
by less than a per cent for variation of the Mi parameter over the range indicated in Fig. |9|. 



4. Conclusions 

Supersymmetric GUTs based on the SO (10) gauge group are highy motivated theories for 
physics beyond the standard model. Many of these models also predict unification of the 
t — b — t Yukawa couplings as well. The restriction of third generation Yukawa coupling 
unification places stringent constraints on the parameter space of supersymmetric models, 
and leads to some new and interesting predictions for new physics. In this paper, we 
restrict ourselves to models with a positive superpotential \i parameter, which is favored 
by BF(b —* S7) and {g — 2)^ measurements. In these models, Yukawa unified sparticle 
mass solutions with R < 1.05 can be found, but only in the HS case, and for very large 
mi6 ~ 5 — 15 TeV case. This also leads generally to TeV scale values of tua and /_t, and 
these numbers, along with TeV scale masses for matter scalars, act to suppress neutralino 
annihilation in the early universe, and hence to predictions of a value of O,^ h ^ ar m excess 
of WMAP allowed values. 
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In this paper, we have explored two methods in which to reconcile Yukawa unified 
sparticle mass solutions with the neutralino relic density. The first case is to allow splitting 
of the third generation of scalars from the first two generations. By decreasing the first and 
second generation scalar masses, we obtain sparticle mass solutions with very light ur and 
cr squark masses, in the 90-120 GeV range. The origin of the light squarks comes from 
the S term in the one-loop RGEs, which is non-zero and large in the case of multi-TeV 
valued split Higgs masses. A search for just two light squarks at the Fermilab Tevatron 
collider using the new Run 2 data should be able to either verify or disprove this scenario. 
In addition, the light squarks give rise to large rates for direct detection of dark matter, 
and should give observable rates for searches at CDMS II. 

Some possible criticisms of the model presented in Sec. |2| is that 1. a degree of fine- 
tuning is required in order to obtain the correct relic density, and 2. the large trilinear 
soft breaking A parameters may lead to charge and/or color breaking minima in the scalar 
potential. We would concur with the first of these criticisms. However, our goal in this 
paper was to examine whether or not a relic density in accord with WMAP could be 
generated in Yukawa unified models with \i > 0, irregardless of fine-tuning arguments. 
Regarding point 2, we note that beginning with Frere et al. [j56|l , bounds have been derived 
on soft breaking parameters of the form \A U \ 2 < 3(m.g + mfj + rn 2 Hu ) by trying to avoid 
charge and/or color breaking minima of the scalar potential. A complete analysis must 
go beyond the tree level and include at least 1-loop radiative corrections to the scalar 
potential |)7]]. We do not include such an analysis here in part because of the philosophical 
question of whether false vacua are allowed, and also because a simple solution of the 
problem may arise by taking the first two generation Aq parameters to be light, and non- 
degenerate with the third generation A parameter. 

The other scenario we examined was the case of non-degenerate gaugino masses at 
the GUT scale. Beginning with any of the Yukawa unified solutions found in an earlier 
study, we find that by dialing M\ to large enough values, a (partially) wino-like LSP can be 
generated, with a relic density in accord with WMAP allowed values. Also, we showed that 
Mi can be dialed so that the value of 2/71^ sits near the Z or h resonance, thus enhancing 
annihilation rates of the Z\ in the early universe. Depending on the model considered, 
this may or may not be enough to bring the neutralino relic density into accord with the 
WMAP measured value. 
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Figure 3: Scan of S0(1Q) HS models with generational non-universality for various planes of 
the parameter space. The colors of the dots correspond to the x 2 value computed from h 2 , 
BF(b -> s-y) and a^ USY . 
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Figure 4: Scan of S'O(IO) models with generational non-universality, in the R vs. fl^ h 2 plane. 
The colors of the dots correspond to the x 2 value computed from h 2 , BF(b — > 57) and a^ USY . 
Points excluded by Edelweiss (CDMS II) direct dark matter searches are denoted by purple (blue) 
dots. 



- 18 - 



1.25 



1.2 



1.15 



1.1 



1.05 




t ? * "% o°°°° 




50 



100 150 200 250 300 350 



m 



u R (GeV) 



Figure 5: Scan of SO(10) HS models with generational non-universality, in the R vs. m^ R plane. 
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Figure 6: Plot of cross section for > 3 jet+ $t events vs. m%, R at the Fermilab Tevatron collider 
after implementing CDF cuts, for same case as in Table [j], but with varying mi6(l). We show 
plots assuming all sparticles contributing to the cross section, and also the cross section coming 
just from gluino and squark production. We also show the 95% CL signal level for errors added 
linearly (upper) or in quadrature (lower). 



- 20 - 



m 



40 



Z^GeV) 

50 60 70 




80 



350 



300 



250 



> 
0) 

200 U 



*3 

150 « 



PS 



100 



50 



Figure 7: Scan of 50(10) HS models with generational non-universality, in the rriu B vs. ra^^ 
plane. The colors of the dots correspond to the % 2 value as listed in Fig. H. 



- 21 - 




Figure 8: Scan of SO(10) HS models with generational non-universality, in the a(Zip) vs. f^/i 2 
plane. The colors of the dots correspond to the x 2 value computed from S7g h 2 , BF(b — > 57) and 
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Figure 9 

Table 1 of Ref. 
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We also show as a horizontal line the relic density limit demanded by WMAP 
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